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Motivated by the LHCb measurement, we analyze the B s -A n + ir~£ + £~ decay in the kinematics 
region where the pion pairs have invariant mass in the range 0.5-1.3 GeV and muon pairs do not 
originate from a resonance. The scalar ■k + tt~ form factor induced by the strange ss current is 
predicted by the unitarized approach rooted in the chiral perturbation theory. Using the two- 
hadron light-cone distribution amplitude, we then can derive the B s -A ty + tt ~ transition form factor 
in the light-cone sum rules approach. Merging these quantities, we present our results for differential 
decay width which can generally agree with the experimental data. More accurate measurements 
at the LHC and KERB in future are helpful to validate our formalism and determine the inputs in 
this approach. 

PACS numbers: 12.39.Fe; 13.20.He; 


Very recently, the LHCb has performed an analysis of rare B s decays into the 7r + 7r /i + /r final state p] and the 
branching fraction is measured as 


B{B S -i = (8.6 ± 1.5 ± 0.7 ± 0.7) x 10“ 


(1) 


where the first two errors are statistical, and systematic respectively. The third error is due to uncertainties on the 
normalization, i.e. the branching fraction of the B° —> J/ip{ —l T + T~)K* (—>■ I\ + ir~). The branching fraction for 
B s -A / o (980W Cl is determined as: 


B{B S -A / O (980)(—17r + 7r-)/xV-) = (8.3 ± 1.7) x 10~ 8 , 


( 2 ) 


which lies in the vicinity of the total branching fraction in Eq. ([!]). Despite the errors, the closeness of the two 
branching fractions and the differential distribution as shown later in Fig. (J4p) may indicate the dominance of the 
/o(980) contributions in the B s -A 7r + 7r“ p + . 

The B s -a- n + n ~is a four-body process. Its decay amplitude shows two distinctive features. On the one side, 
the 7r + 7r - final state interaction is constrained by unitarity and analyticity. On the other side, the b mass scale is much 
higher than the hadronic scale Aqcd, which allows an expansion of the hard-scattering kernels in terms of the strong 
coupling constant a s and the dimensionless power-scaling parameter Aqcd /nib- In Refs. PHI], we have developed 
a formalism that makes use of these two advantages. This approach was also pioneered in Ref. B10. and see also 
Refs. HHII1 for applications to charmless three-body B decays. In doing this, the new formalism can simultaneously 
merge the perturbation theory at the mj, scale and the low-energy effective theory based on the chiral symmetry to 
describe the S-wave mr scattering. The aim of this work is to further examine this formalism by confronting this 
theoretical framework with the recent data on B s -A 7r + 7r _ /i + /i - . An independent analysis that is based on the 
perturbative QCD approach is also under progress )12j 

We start with the differential decay width for B s -A n + Tr~£ + £~. The effective Hamiltonian for the transition 
b -A s£+£~ 
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involves various four-quark and the magnetic penguin operators Oj. the Ci(/i) are the corresponding Wilson coefficients 
for these local operators Oi. Gf is the Fermi constant, and Vtb = 0.99914 ± 0.00005 and V ts = —0.0405to;oi2 ECS are 
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the CKM matrix elements. The b and s quark masses are m b = (4.66 ± 0.03)GeV and m s = (0.095 ± 0.005)GeV |13| . 
The b —> sl + l~ transition has the decay amplitude 


iM(b -»• stl~) = iN r x j(C 9 + C w )[sb\ V - A [M]v+A + (C 9 - C 10 )[sb} V -A[U}v-A 

+ 4 C 7 L m b [sia lx „(l + 7s)fr]^y x [^7^] + 4C , 7HTO 6 [sicr jUy (l -75)6]^- x [£7^] }>, 


where CVl = C 7 and C 7 r = C 7 Lm s /mb, and 


Ni = ^°^v tb v; s . 


(3) 


(4) 


4^2 7T 

The B —► MiM 2 l + £~ is a four-body decay mode, whose decay amplitude can be obtained by sandwiching Eq. ([ 3 ]) 
between the initial and final hadronic states. The spinor product [sb] will be replaced by corresponding hadronic 
matrix elements. A general differential decay width for B —> M-\ M 2 t + t~ with various partial wave contributions has 
been derived using the helicity amplitude in Ref. [14] . In the B s —> p + pT case, the S-wave contribution will 

dominate and thus the angular distribution is derived as 

d 3 T 3 r 


= g[JT +J 2 c cos( 20 O 


(5) 


dm 2 v dq 2 d cos 9i 

where 9i is the polar angle between the p~ and the B s moving direction in the lepton pair rest frame. The angular 
coefficients are given by 

( 6 ) 
(7) 


Ji = |l-4ior + l-4/?or + 8mf |^2, 0 Ar 0 | cos (^lo ~ $ro) + 4f?q |Aq 
J 2 = ^A 2 ||-4lo| 2 + l-T/sol 
In the above equations, /3^ = \Jl — 4 m 2 /q 2 , and me = me/\fq 2 ■ The helicity amplitude is 
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where 




( 8 ) 


(9) 


Here the script t denotes the time-like component of a virtual state decays into a lepton pair. The function A is 
related to the magnitude of the 7r + 7r _ momentum in B s meson rest frame: A = A (m 2 B , 77i 2 +7r _ , q 2 ), and A (a 2 , fe 2 , c 2 ) = 
(a 2 — 6 2 — c 2 ) 2 — 46 2 c 2 . The combination of the time-like decay amplitude is introduced in the differential distribution 


A/ = A% t - A° Lit = 2 y/W 2 C, 


10*- 
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( 10 ) 


7T7T form factors used in Eq. ([8]) are defined by 
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The B s —» 
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(( tt+tt )s(Pnir)\sa fJ ,„q , '"/ 5 b\B s (pB s )) = 


FAml^.q 2 ) 


(m 2 B - m 2 Jg M - q 2 Pfi 


( 11 ) 


m nn (m B +m 

As we have shown in Ref. [2], an explicit calculation of the B s —> n + n~ form factors requests the knowledge on 
generalised light-cone distribution amplitudes [T6H22]- The expressions in the light-cone sum rules are given as [2j, 
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TABLE I: The B s —> /o(980) form factors in the light-cone sum rules at LO and NLO in a a [15] . 


LO F(0) a F b F 

NLO F(0) a F b F 

Fi 0.185 ±0.029 1.44±g;J| 0.59ig;^ 
F 0 0.185 ± 0.029 0.47t°;Jg O.Oltg.dd 
Ft 0.228 ±0.036 1.42±g;^g 0.60^;^ 

T\ 0.238 ±0.036 1.50±g;Jl 0.58ig;S? 

F 0 0.238 ± 0.036 0.53ld;io —0.36±o Sl 
Ft 0.308 ±0.049 1.46t(U °- 58 ±o.o 7 
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where 


»r D p / 2 xTO b + m s 

tv F = B 0 F^{m^) 2 m 2 j b ex P 
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the above the scalar irn form factor is defined as 


L ^ J 

< 2 - So + V(min + 9 2 - gp) 2 + 4m^(mg - g 2 ) 
2m; 2 


( 16 ) 


( 0 | ss | f + f ) = B 0 F^(ml n ), ( 17 ) 

and the Bg is the QCD condensate parameter: 

<0 |ot|0 ) = -/ 2 Bo, (18) 

with MeV being the pion decay constant at LO. For the numerics, we use f n = 91.4MeV and (0|qq|0) = —(0.24± 
0.01)GeV 3 (for a review see Ref. [23]), which corresponds to Bq = (1.7 ± 0.2) GeV. The M is a Borel parameter 
introduced to suppress higher twist contributions. Our formulae can be compared to the results for the B s —► / 0 (980) 
transition HE with the correspondence 


to/ 0 <± m„, 


K 0 (n) ++ KM, ff 0 ++ B 0 F^(ml,), 


( 19 ) 
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FIG. 1: Feynman diagrams for the scalar form factor at tree-level and one-loop level in CHPT. The wave function renormal¬ 
ization diagrams are not shown here. 


where ff 0 is the decay constant of / 0 (980) defined by the scalar current. The twist-3 distribution amplitudes, 
and 4for the scalar 7T7 t state have the same asymptotic forms with the ones for a scalar resonance [23], while 
the twist ones can be similarly expanded in terms of the Gegenbauer moments. Inspired by this similarity, we can 
plausibly introduce an intuitive matching: 

Jf q 2 ) = ^B 0 F^{ml n )F^ fa {q 2 ). (20) 

J fo 

Here we have assumed the dominance of the /o(980) which is justified in the B s —> 7r + 7r _ /i + /i _ as shown in the data 
in Eq. ([2]) and Eq. Q. 

The B s —> / o (980) form factors have been calculated in the light-cone sum rules at leading order (LO) and next-to- 
leading order (NLO) in ot s [T51 I25H2T] . and in the perturbative QCD approach 12M55] in Ref. [MJ. The momentum 
distribution in the form factors has been parametrized in the form: 


Fitf) 


_ pm _. 

1 - a.iq 2 /m B f + bi{q 2 /m 2 Bs )' 


( 21 ) 


Numerical results for these quantities where // 0 = (0.18 ± 0.015) GeV [35] are taken from Ref. [15] and are collected 
in Tab.JIj Using a different value for ff 0 for instance in Ref. HUES] will not induce any difference to the generalized 
form factor, since such effects will cancel as demonstrated in Eq. (20). In the following calculation, we will use the 
NLO results for the B s —> f 0 transition. Using the LO results can reduce the differential decay width by about 40%. 

The scalar tttt form factor F nw (ml n ) has been calculated within a variety of approaches using (unitarized) chiral 
perturbation theory (CHPT) [37H44] and dispersion relations [45]. In terms of the isoscalar S-wa,ve states 


= ^fl 7r+7r ) + ^ I^V 0 ), 

\KK) 1=0 = | K+K~) + ±= | K°K°) , 

the scalar form factors are defined as 

V2B 0 Ff(s) = (0|ss|7T7r)i =o , 
V2B 0 F$(s) = {0\ss\KK) l=o , 


( 22 ) 

(23) 


(24) 

(25) 


where the notation (tt = 1, K = 2) has been introduced for simplicity, and the convention F n7r (ml n ) = 2/a/3 Fl(m 2 % ). 
In the CHPT, expressions have already been derived by calculating the diagrams in Fig. [I] up to NLO 

16j?j2 So o 0 ™ 2 


j7’CHPT( g ) = ^ 


[ P 


8 c e O 777 ‘ 

(2 LI - LI) + ~p.L\ + ^J r KK {s) + g -fW 
4 LI 


F 2 chpt (s) = 1 + yr( s - m l- 4 m 2 K ) + yi(s- 4m 2 K ) + (4m 2 K + ml) + 


32 LI 


p ' m K + gAO 


9s — 8 m 2 K 
18 P 


p v i s ) + Jkk( s )- 


(26) 


(27) 
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FIG. 2: The s-channel diagrams to the scalar nir form factors in CHPT. With the increase of the 7T7 t invariant mass, higher-order 
contributions may become important. In the unitarized approach [46], these diagrams can be summed. 


With the increase of the invariant mass of the irn system, higher order contributions become more important. It 
has been proposed that the unitarized approach which can sum higher order corrections and extend the applicability 
to the scale around 1 GeV [32]. A sketch of the resummation scheme is shown in Fig. [2j In this figure, the K(s) is 
the A-wave projected kernel of meson-meson scattering amplitudes [421 1431 : 


K(s) 


Ku I< 12 

K 2 i K 22 


K n = 


2s - nil 

2 f 2 ’ 


Kvi = K-n 


\/3s 3s 

4J2 ’ A22 = 4^2 > 


(28) 

(29) 


where the subscript 1,2 denotes the 7T7 t and KK state respectively. The function g{s) is the loop integral which can 
be calculated in the cutoff-regularization scheme with q max ~ lGeV being the cutoff [cf. Erratum of Ref. [46]] or in 
dimensional regularization. In the latter scheme, the meson loop function gu(s) is given by 


Jliis) 


1 

167T 2 


1 - log 



CTi(s) log 


f o-i(s) + i y 

Vcr^s) - 1 )_ 


-gu(s). 


(30) 


with (Tj(s) = — 4 ml/s. Imposing the unitarity constraints, the scalar form factor can be expressed in terms of the 

algebraic coupled-channel equation [SB] US] 


F(s) = R(s)[I + jtsjfffs)]- 1 

= i?(s)[/- 5 (s)A7s)] + 0(p 6 ), 


(31) 


where the above equation has been expanded up to NLO in the chiral expansion. The R(s) = (R±(s), R 2 (s)) includes 
both tree-level contributions, and other higher order corrections that have not been summed. Thus this function has 
no right-hand cut, and can be obtained by matching onto the CHPT results in Eqs. (26|27l [301 Wf\ \ 


Ri{s) 

Ms) 


\/3 j 16 ml 


7,2 Sc 

8 LI 


72tt 2 /2 




AT r ^f\T r 

1 + (s - 4 m% -mD + ^pis^ 4m|) + ^ (4 m* K + ml) 


L K 


36tt 2 /2 


1 T log 



32 LI 2 2 

P m K + 2 ^ 1 ) 


(32) 


(33) 


With the above formulae and the fitted results for the low-energy constants L [ in Ref. [JO] (evolved from m p to the 
scale /i = 2g max /y / e), we show the strange 7T7 t form factor in Fig. [3j The modulus, real part and imaginary part are 
shown as solid, dashed and dotted curves. 

Equipped with the results for scalar form factor and heavy to light transition, we can explore the differential 
branching fraction for the B s —> ir + n~ g + . Our theoretical results for dB/dm^r is given in the left panel of Fig. [4j 
This clearly shows the peak corresponding to the /o(980). In order to compare with the experimental data, we also 
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FIG. 3: The nn scalar form factor obtained in the unitarized chiral perturbation theory. The modulus, real part and imaginary 
part are shown in solid, dashed and dotted curves. 
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FIG. 4: The differential branching ratio for the B s —> n + n~£ + £~. The experimental data (with triangle markers) has been 
normalized to the central value of the branching fraction: B(Bg tt + 7t - p + p“) = (8.6 ± 1.5 ± 0.7 ± 0.7) x 10 _s . Theoretical 
predictions (with square markers) are based on the result for the time-like scalar form factors derived in the unitarized CHPT. 


give the binned results on the right panel in Fig. [4] from 0.5 GeV to 1.3 GeV. Theoretical errors shown in this panel 
arise from the ones in the form factors. The experimental data (with triangle markers) has been normalized to the 
central value given in Eq. 0. The comparison in this panel shows a general agreement between our theoretical 
prediction and the experimental data except in a few bins. This agreement is very encouraging. 

In spite of the agreement, there exist some differences in our results and data. For instance our theoretical result 
does not show the enhancement at m nn ~ (800,1100,1250)MeV as given in the data. The excess at 800 MeV may 
come from the tail of the B s —> ij (—» 7T + 7r”7r°, ir + ir~”/)/j, + , while in the range above lGeV, the contribution from 
the /o(1370) may not be negligible. 

Integrating out the m we have the branching fraction: 

B(B S -> / O (980)(—>■ tt + tt-)h + h~) = (4.1 ± 1.6) x 10~ 8 , (34) 

which deviates from the data by about 2er. However, one expects the experimental result in Eq. ([2]) would get 
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somewhat reduced. This can be witnessed by the B —> J/'ipK and B —tK n + /j, [ T3] 

(4.49 ± 0.23) x 10" 7 _ 4 ^ 

B(B~ —>• J/4’K~) (1.027 ± 0.031) x lO- 3 ' { ! 

If this ratio were not sensitive the light meson in the final state which is true in most cases, the branching fraction 
for the B s —1 J/ipf o (980) [13] 


B(B S -*• J/tpfo) = (1.39 ± 0.14) x 1(T 4 , 


would indicate 


B(B S -)■ /o(980)/r + /r _ ) - 6.1 x 1CT 8 . (36) 

This value is smaller by about 30% than the central value given in Eq. 0 , and is more consistent with our theoretical 
result. The future measurement with more data at the experimental facilities like LHC and KEKB will be able 
to clarify this point, and thus to examine our theoretical formalism more precisely. We strongly encourage our 
experimental colleagues to conduct such measurements. 

In summary, in this work we have analyzed the B s —> tt + tt~£ + £~ that has focused on the region where the pion pairs 
have invariant mass in the range 0.5-1.3 GeV and muon pairs do not originate from a resonance. We have adopted the 
approach proposed in our previous work EHS (see also Ref. [5j for an overview) which makes uses of the two-hadron 
light-cone distribution amplitude. The scalar 7r + 7r _ form factor induced by the strange ss current is predicted by 
the unitarized chiral perturbation theory. The heavy to light transition can then be handled by the light-cone sum 
rules approach. Merging these quantities, we have presented our theoretical results for differential decay width and 
compared with the experimental data. Except in a few bins, our theoretical results are in alignment with the data. 
We have also discussed the disagreement and given our expectation. More accurate measurements at the LHC and 
KEKB in future are helpful to validate/falsify our formalism and determine the inputs in this approach. 
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